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Abstract 

Background: Increased airway wall thickness {AWT) and parenchymal lung destruction both contribute to airflow 
limitation. Advances in computed tomography (CT) post-processing imaging allow to quantify these features. The 
aim of this Dutch population study is to assess the relationships between AWT, lung function, emphysema and 
respiratory symptoms. 

Methods: AWT and emphysema were assessed by low-dose CT in 500 male heavy smokers, randomly selected 
from a lung cancer screening population. AWT was measured in each lung lobe in cross-sectionally reformatted 
images with an automated imaging program at locations with an internal diameter of 3.5 mm, and validated in 
smaller cohorts of patients. The 15 th percentile method (Perc15) was used to assess the severity of emphysema. 
Information about respiratory symptoms and smoking behavior was collected by questionnaires and lung function 
by spirometry. 

Results: Median AV\AT in airways with an internal diameter of 3.5 mm (AWT 35 ) was 0.57 (0.44 - 0.74) mm. Median 
AWT in subjects without symptoms was 0.52 (0.41-0.66) and in those with dyspnea and/or wheezing 0.65 (0.52-0.81) 
mm (p<0.001). In the multivariate analysis only AWT 35 and emphysema independently explained 31 .1 %and 9.5%of 
the variance in FEV^/opredicted, respectively, after adjustment for smoking behavior. 

Conclusions: Post processing standardization of airway wall measurements provides a reliable and useful method 
to assess airway wall thickness. Increased airway wall thickness contributes more to airflow limitation than 
emphysema in a smoking male population even after adjustment for smoking behavior. 
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Introduction 

The quantification of airway dimensions by CT has be- 
come feasible with the development of multi detector 
computed tomography (CT) and new software tools for 
image analysis [1,2], Assessment of airway dimensions 
by CT has been studied particularly in relation to 
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asthma, smoking and chronic obstructive pulmonary 
disease (COPD) [3-8], diseases generally associated with 
chronic or intermittent airflow limitation. So far, airway 
wall thickness (AWT) measurements have been per- 
formed by selecting well quantifiable airways [9-12] or 
by standardizing dimensions to airways with a 10 mm 
internal lumen perimeter (pi 10, equivalent to about 
3.2 mm internal airway lumen diameter) derived from a 
small number of airways [5,6,13-15]. 

More recently, low-dose multi slice CT has become 
available, a technique with a good quantifying performance 
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that is preferred for monitoring of pulmonary and airway 
abnormalities as compared to the high radiation exposures 
with high-resolution computed tomography (high-dose 
CT). The cumulative radiation dose exposure with low- 
dose CT remains very low, even when individuals are 
exposed multiple times [16]. Airway dimensions measured 
with low-dose CT have only been reported in few studies 
using the same diversity in analytic approaches as applied 
with high-dose CT measurements [7,17-19]. 

Airway dimensions and extent of emphysema are known 
to be associated with airflow limitation [6,9-14,19,20], 
although the influence of smoking behavior or signs of 
airway disease such as cough, dyspnea, wheezing or 
mucus overproduction on airflow limitation is not clear 
[5,7,21,22]. 

The aims of this study are to quantify airway dimen- 
sions of the lung in multiple airway sections of each lobe 
in a novel manner and the extent of emphysema by 
using low-dose CT. These measurements were asso- 
ciated with the influence of airflow limitation, respira- 
tory symptoms and corrected for smoking behavior. 

Methods 

Population 

We randomly selected 500 current and former smokers 
participating in the Groningen cohort of a male 
population-based multi-centre lung cancer screening study 
(NELSON). The Dutch ministry of health and the Medical 
Ethics Committee of the hospital approved the study 
protocol. Informed consent was obtained from all partici- 
pants. Detailed inclusion criteria and characteristics have 
been described elsewhere [23]. In short, subjects with a 
smoking history of at least 20 pack-years were included. 
Information on the presence or absence of respiratory 
symptoms and smoking (pack-years and former or current 
smoking) was obtained by questionnaires. The question 
used to record respiratory symptoms was "do you have 
experienced the following symptoms cough, sputum 
expectoration, wheezing or dyspnea for at least 3 months 
during the past year, even when you did not have a cold?" 

Lung function 

All participants performed standardized spirometry accor- 
ding to the European Respiratory Society guidelines [24], 
including forced expiratory volume in 1 sec (FEVi) and 
forced vital capacity (FVC) at the start of the study. In 
this population- based study we did not administer a 
bronchodilator. 

CT scanning 

Low- dose CT images of the chest were acquired at full 
inspiration after appropriate instruction on one CT scan- 
ner (Sensation 16, Siemens Medical Solutions, For- 
chheim, Germany) [23] according to the following 



protocol: spiral acquisition at 120 kV, 20mAs, rotation 
time 0.5 s, pitch 1.5 and collimation 16x0.6 mm, field of 
view 300 to 350 mm, slice thickness 1 mm and slice 
increment 0.7 mm. The effective radiation dose was less 
than 0.8 mSv. Contrast medium was not used. The 
images were reconstructed to a pixel matrix of 512x512 
using B30f kernel. Thus, the spatial resolution was 0.59 
to 0.68 mm in x-y plane, and 0.7 mm in z plane. The CT 
system was calibrated routinely. 

Quantification of AWT 

AWT was measured in cross-sectionally reformatted 
images with an automated research software prototype 
MEVIS Airway Examiner vl.O (release 2009, Fraunhofer 
MEVIS, Bremen, Germany) based on an algorithm by 
Weinheimer at locations with a fixed internal diameter 
of 3.5 mm in each lung lobe [25]. This software auto- 
matically extracts airway centerlines, re-samples images 
perpendicular to the airway direction at equally spaced 
positions along the centerline and detects inner and 
outer airway wall borders in these images. The outer 
wall border is detectable when no adjacent tissue with 
similar CT density is present and is taken into account 
when the wall is detected in at least 25% of the perim- 
eter at a location. AWT and the fraction of perimeter 
where the outer wall border was identified (Assessed 
Perimeter Fraction, APF) are calculated for each loca- 
tion. Wall thickness quantification accounts for partial 
volume effects by integrating Hounsfield units across the 
wall. Accuracy and reproducibility of this algorithm 
was tested previously under clinical conditions using a 
similar protocol as used in our study [2]. Average wall 
thickness and cumulative APF of all detectable airway 
locations with a fixed lumen diameter is reported per 
lobe and for the whole lung. The borders of the lung 
lobes were automatically calculated by the software in 
a standard way. All low-dose CT scans were visually 
evaluated for appropriate segmentation. 

Quantification of emphysema and lung volume 

Quantification of emphysema was based on density diffe- 
rences and measured with a software tool called Image 
Xplorer (Image Sciences Institute, Utrecht, the Nether- 
lands) [16,26]. This software produces automatically the 
lung volume. The extent of emphysema was automatically 
performed at the 15 th percentile (Percl5) of the Houns- 
field density distribution. Percl5 is the threshold density 
value where 15% of all voxels has a lower density [27]. A 
lower Percl5, i.e. closer to -1000 HU, means that more 
emphysema is present. All scans were reconstructed with 
a soft reconstruction filter (Philips B, Siemens B30f). Air- 
ways were automatically excluded to assess density of lung 
parenchyma exclusively and HU densities of the entire scan 
were recalibrated using automatically measured average 
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densities in the trachea and shifting the HU values of the 
entire scan so that air density in the trachea became -1000 
HU. Additionally, the percentage of low attenuation area, 
defined as the proportion of low-density voxels below -950 
HU (%LAA-950HU) was used. %LAA-950HU was log- 
transformed because of skewed distribution. 

Explorative studies 

Prior to the research described above we have 

1) established the optimal internal airway diameter, i.e. 
the internal airway diameter that allows the highest 
number of cumulatively Assessed Perimeter 
Fractions (APF) for the whole lung. Therefore we 
measured APF on 20 selected NELSON CTs in 
airways with a lumen diameter of 2.5, 3.0, 3.5, 4.0, 
4.5 and 5.0 mm (± 0.25 mm) divided into 3 groups: 
no emphysema and normal lung function (n = 8, p 15 

> -920 and FEVJFVC > 85 %), moderate 
emphysema and normal lung function (n = 4, -940 < 
p 15 < -960 and FEV^FVC > 70 %) and no 
emphysema and severe airflow limitation (n = 8, p 15 

> -920 and FEVJFVC < 50 %). 

2) compared the mean AWT 3 5 , using the same method 
as described above, at the established internal lumen 
size with high- and low-dose CT in 8 NELSON 
subjects from whom high- and low-dose CT were 
available. These CT data were obtained in spiral 
mode with 16x0.75 mm collimation and in full 
inspiration with the same scanner (Sensation- 16 
Siemens Medical Solutions, Forchheim, Germany). 
Axial images were reconstructed with 1.0 mm 
thickness at 0.7 mm increments. All scans were 
reconstructed with a soft reconstruction filter 
(Siemens B30f) at a 512x512 matrix. 

3) determined the generation where airways with the 
established optimal internal lumen size are present. 
AWT measurements at 3.5 mm internal lumen size 
were performed in 57 randomly selected low-dose 
CTs of NELSON subjects. A multi-planar 
reconstruction (MPR) was made in each case of the 
apical upper lobe bronchus (Bl) and the posterior 
lower lobe bronchus (B10). Subsequently the 
location was projected on the segmentation image. 
Three-dimensional image moving created the 
opportunity to observe airways from various 
directions and to check bifurcations and count 
airway generations according to the method of 
Boyden [28]. 

Statistical analysis 

Data are reported as mean ± standard deviation (SD) or 
median (25 th - 75 th percentile) values as appropriate. 
The mean AWT at 3.5 mm internal lumen size 



(AWT3.5) of all five lobes per case was calculated taking 
into account the APF per lobe by the following formu- 
lae: ((AWT left upper lobe x APF left upper lobe) + 
(AWT left lower lobe x APF left lower lobe) + (AWT 
right upper lobe x APF right upper lobe) + (AWT right 
middle lobe x APF right middle lobe) + (AWT left 
upper lobe x APF left upper lobe)) / sum of APF of all 
lung lobes. AWT 3 5 for the whole population was skewed 
distributed, therefore we report median AWT and range, 
and log-transformed AWT was used in the analyses. 

Differences between various categories were explored 
using chi-square tests (dichotomous data), 2-tailed 
unpaired Students t-tests for normally distributed con- 
tinuous data and Mann- Whitney U-tests for not nor- 
mally distributed continuous data. The difference in 
airway wall thickness between lung lobes was assessed 
with a Wilcoxon signed rank test. Univariate linear re- 
gression analyses was used to study the relationships 
between clinical variables and AWT, and those variables 
with FEV!%predicted. Next, multivariate linear regres- 
sion analyses were performed on those clinical variables 
showing significance in the univariate regression ana- 
lyses. Outcomes of these analyses have been described 
with betas and p-values. Bland- Altman plot was used to 
analyze the agreement between AWT by high- and low- 
dose CT [29]. All statistical analyses were performed 
using SPSS 18.0 for Windows; P-values below 0.05 were 
considered statistically significant. 

Results 

Population characteristics 

After visual evaluation 8 out of the 500 randomly selected 
subjects were excluded because of (partial) missing of air- 
way segmentation on CT. The mean age of the cohort was 
59.4 (± 5.2) years, approximately 59% were current smo- 
kers and median pack-years was 34.0 (28.0 - 45.6). More 
than 51% of the population reported at least one respira- 
tory symptom (Table 1). 

Airway wall thickness 
Exploratory analyses 

The highest numbers of cumulatively assessed perimeter 
fractions (APF) of airways were reached at an internal 
lumen perimeter of 3.5 ± 0.25 mm (Figure 1); therefore 
this diameter was selected for further analysis. Median 
AWT3.5 on low-dose CT was comparable with median 
AWT3.5 on high-dose CT, respectively 0.57 (0.48 - 0.74) 
mm and 0.55 (0.47 - 0.73) mm (p = 0.89, n = 8). This 
demonstrates that MEVIS software analysis of data from 
low-dose CT gives similar results as from high-dose CT 
(Figure 2). Airways with an internal diameter of 3.5 mm 
appeared mainly in the 5 th generation (range, 3 rd - 7 th 
generation) in the upper lobe bronchus, and mainly in 
the 8 th - 10 th generation (range, 6 th - 12 th generation) in 
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Table 1 Clinical and demographic characteristics of heavy 
smokers from the general population cohort 

N 

Age, years 
Pack-years smoking 
Current smoking,% 
FEV 1; liter 
FEV 1; %predicted 
FEVFVQo/o 

Emphysema; Perc15, HU 
Emphysema;%LAA -950 HU 
Emphysema; >5%LAA -950 HU,% 
Lung volume on CT, liter 
Chronic Mucus Hypersecretion/^ 
Cough,% 
Dyspnea,% 
Wheezing,% 

No respiratory symptoms,% 



492 

59.4 ± 5.2 

34.0 (28.0 - 45.6) 

59.1 

3.45 ± 0.76 

98.2 ± 19.7 

70.0 ± 10.7 

-920 (-930 to -907) 

2.5 (1.3-4.3) 

19.7 

6.5 ± 1 .4 

29.7 

32.7 

28.4 

25.3 

48.6 



Mean ± standard deviation shown for continuous data and median 
(interquartile range) for non-parametric distribution. 
Definition of abbreviations: HU = Hounsfield Units; Perc15 = the threshold 
density value where 15% of all voxels has a lower density; %LAA -950 HU = 
percentage low attenuation areas < -950 HU; >5%LAA -950 HU,% = 
percentage of the population having >5% low attenuation areas < -950 HU. 

the lower lobe bronchus. This distribution was observed 
irrespective of smoking, presence of airflow limitation 
(defined as FEXVFVC <0.70) or chronic mucus hyperse- 
cretion (CMH) (Figure 3). 



no emphysema, normal lung function 
moderate emphysema, normal lung function 
no emphysema, severe airflow limitation 



> < 



E ~ 




3.0 3.5 4.0 4.5 

internal lumen diameter (mm) 

Figure 1 Determination of the optimal airway size. Mean 
cumulative assessed perimeter fractions in the total lung, for 
different groups of patients at different internal airway lumen 
diameters. APF was measured on low-dose CT of 20 selected 
NELSON subjects divided into 3 groups: subjects without 
emphysema and with normal lung function (n = 8, perc15 > -910 
and FEVyFVC > 85 %), with moderate emphysema and normal lung 
function (n = 4, -940 < perc15 < -960 and FEVYFVC > 70 %) and 
without emphysema and having severe airflow limitation (n = 8, 
perc15 > -920 and FEVYFVC < 50 %), in airways with a lumen 
diameter of 2.5, 3.0, 3.5, 4.0, 4.5 and 5.0 mm (± 0.25 mm). 



Airway wall thickness in the population 

In the whole population the median (25 th - 75 th percent- 
ile) AWT3.5 was 0.57 (0.44 - 0.74) mm. The APF in the 
whole lung varied from 142 to 295 (median 215). The 
results per lung lobe are presented in the online 
Additional file 1: Table SI). 

Airway wall thickness and clinical characteristics 

Significantly higher AWT 35 values were observed in 
subjects with dyspnea and/or wheezing (n = 181, median 
AWT3.5 0.66 mm) or with cough and/or CMH (n = 201, 
median AWT 3 5 0.63 mm) compared to subjects without 
dyspnea and/or wheezing (n = 309, median AWT 3 5 
0.53 mm, p<0.001) or without cough and/or CMH 
(n= 291, median AWT 0.53 mm, p-values < 0.001). 
Current smokers and former smokers had comparable 
median AWT 35 values, i.e. 0.58 mm and 0.56 mm 
(p = 0.113) respectively (Table 2). 

Univariate linear regression analysis showed inverse 
associations between log-AWT 3 5 and FEVi (b = -0.233, 
p < 0.001), FEVYFVC (b = -0.015, p < 0.001), FEV^/o 
predicted (b = -0.010, p < 0.001) and lung volume 
(b = -0.055, p<0.001) and positive associations between 
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Figure 2 Comparison of airway wall thickness on high- and 
low-dose CT. AWT 3 . 5 was measured on high and low-dose CT of 8 
NELSON subjects. Bland & Altman plot shows agreement between 
mean AWT 35 determined by high- and low-dose CT. Dashed lines 
depict the 95% confidence interval. 
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Figure 3 Distribution of 3.5 mm sized airways over airway generations. The distribution of 3.5 mm internal lumen sized airways over the 
2 nd till 7 th airway generations in the apical upper lobe bronchus (B1) and over the 6 th till 12 th generation in the lower lobe bronchus (B10) in the 
right lung, for subjects with and without airway obstruction, for subjects with and without Chronic Mucus Hypersecretion (CMH) and for current 
and former smokers. The number of cases (e.g. no CMH and CMH) in each group was limited but the distribution was similar. 



log-AWT 3 5 , Percl5 and number of pack-years smoking, 
respectively (b = 0.003, p < 0.001) and (b = 0.003, 
p = 0.003) (Table 3). 

Multivariate analysis on all clinical variables signifi- 
cantly associated with AWT in univariate analyses 
revealed that log-AWT 35 was independently asso- 
ciated with lower FEV^predicted (b = -0.010, 
p < 0.001), higher Percl5 (b = 0.005, p < 0.001), and 
lung volume (b = -0.037, p = 0.005) respectively 
(Table 4). 



Contribution of airway wall thickness and emphysema to 
airflow limitation 

To study the contribution of AWT 35 , emphysema 
(Percl5), and clinical variables to airflow limitation, 
FEVi%predicted was taken as dependent variable. A sig- 
nificant negative association was found between FEV^ 
predicted and log-AWT 3 5 (b = -31.21, p < 0.001), pack- 
years (b = -0.153, p = 0.004), cough, CMH, dyspnea and 
wheezing (b = -9.234, -8.17, -11.99, -13.26 respectively, 
all p-values < 0.001), and a positive association between 



Table 2 Airway wall thickness in subjects with and without respiratory symptoms and in current and former smokers 





n 


AWT3.5 (mm) 




n 


AWT3.5 (mm) 


p-value 


CMH 


146 


0.62 (0.49-0.80) 


no CMH 


346 


0.56 (0.44-0.71) 


0.002* 


Cough 


161 


0.64 (0.50-0.84) 


no cough 


331 


0.54 (0.43-0.70) 


<0.001* 


Dyspnea 


140 


0.65 (0.51-0.80) 


no dyspnea 


352 


0.54 (0.43-0.69) 


<0.001* 


Wheezing 


125 


0.66 (0.52-0.85) 


no wheezing 


367 


0.53 (0.43-0.70) 


<0.001* 


CMH and/or cough 


201 


0.63 (0.49-0.79) 


no CMH or cough 


291 


0.53 (0.43-0.70) 


<0.001* 


Dyspnea and/or wheezing 


181 


0.66 (0.52-0.83) 


no dyspnea or wheezing 


311 


0.53 (0.42-0.67) 


<0.001* 


Cough, CMH, dyspnea and wheezing 


49 


0.69 (0.51-0.89) 


no cough, CMH, dyspnea or wheezing 


239 


0.52 (0.41-0.66) 


<0.001* 


Current smoking 


291 


0.58 (0.45-0.75) 


Former smoking 


201 


0.56 (0.43-0.71) 


0.113 



Data presented as median (interquartile range) values. 

Definition of abbreviations: AWT 3 . 5 = airway wall thickness at 3.5 mm internal lumen size; CMH = chronic mucus hypersecretion. 
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Table 3 Univariate associations between (A) log-transformed AWT 3 . 5 and (B) FEV n % predicted, and clinical 
characteristics 



Dependent variable A. Log-AWT 3 5 B. FEV l7 % predicted 





Beta 


p-value 


Beta 


p-value 


Log-AV\^3. 5 






-31.21 


<0.001 


FEVt % predicted 


-0.010 


<0.001 






FEV 1; liter 


-0.233 


<0.001 


23.56 


<0.001 


FEVFVQo/o 


-0.015 


<0.001 


1.441 


<0.001 


FVC,% predicted 


-0.009 


<0.001 


0.971 


<0.001 


Emphysema; Perc15 


0.003 


<0.001 


0.204 


<0.001 


Emphysema; Log-%LAA -950 HU 


-0.041 


0.023 


-5.872 


<0.001 


Lung volume on CT, Liter 


-0.055 


<0.001 


-0.257 


0.691 


Pack-years 


0.003 


0.003 


-0.153 


0.004 


Smoking (former/current) 


0.059 


0.065 


-3.257 


0.071 


Age, years 


-0.002 


0.458 


-0.289 


0.092 


Height, cm 


-0.002 


0.444 


-0.121 


0.425 


Cough 


0.156 


<0.001 


-9.234 


<0.001 


CMH 


0.114 


0.001 


-8.171 


<0.001 


Dyspnea 


0.158 


<0.001 


-11.99 


<0.001 


Wheezing 


0.202 


<0.001 


-13.26 


<0.001 



Definition of abbreviations: Log-AWT 3 5 = log transformed airway wall thickness at 3.5 mm diameter; CMH = chronic mucus hypersecretion; Log-%LAA -950 
HU = log transformed percentage of low attenuation areas < -950 HU. 



FEV%predicted and Percl5 (b = 0.204, p < 0.001) 
(Table 3). 

Multivariate analysis on all clinical variables that were 
associated with FEV^predicted in univariate analyses 
showed that higher FEV^predicted was significantly 
associated with lower AWT 35 values (b = -31.3, 
p < 0.001) and with higher Percl5 (b = 0.342, p < 0.001) 
(Table 4). AWT 3 . 5 and Percl5 explained 31.1% and 9.5% 
of the variance of FEV^predicted, respectively. The 



Table 4 Multivariate linear regression: dependent 
variable is (A) log transformed AWT 3 . 5 and (B) FEV^/o 
predicted 



Dependent variable 


A. Log-AWT 3 . 5 


B. FEV l7 % 


predicted 


Beta 


p-value 


Beta 


p-value 


FEV 1; % predicted 


-0.010 


<0.001 






Log-AV\^3. 5 






-31.28 


<0.001 


Emphysema; Perc15 


0.005 


<0.001 


0.342 


<0.001 


Lung volume 


-0.037 


0.005 


0.209 


0.773 


Pack-years 


0.002 


0.031 


-0.029 


0.499 


Smoking (former/current) 


-0.01 1 


0.700 


-2.756 


0.086 


Cough 


0.030 


0.404 


-0.408 


0.834 


CMH 


-0.026 


0.447 


-2.098 


0.262 


Dyspnea 


-0.007 


0.828 


-2.945 


0.116 


Wheezing 


0.058 


0.124 


-4.004 


0.051 



Definition of abbreviations: Log-AWT 3 5 = log transformed airway wall thickness 
at 3.5 mm diameter; CMH = chronic mucus hypersecretion. 



results of the multivariate regression analysis with em- 
physema expressed as %LAA-950 HU as independent 
variable, are presented in the online Additional file 1: 
Table S2. 

Discussion 

Low-dose CT is an appealing approach to quantify si- 
multaneously pulmonary and airway abnormalities. Our 
study shows that the use of low-dose CT combined with 
modern post processing software tools provides reliable 
information on airway wall thickness and the extent of 
emphysema in a heavy smoking male population. Al- 
though CT does not provide dynamic measurements, 
airway wall thickening and emphysema explained re- 
spectively 31.1% and 9.5% of the variance in FEV^pre- 
dicted, the most commonly used variable of airflow 
limitation. Changes in AWT of more than 0.1 mm 
reflecting lumen surface changes over 8% measured at 
one air lumen level were observed between cases with 
and without respiratory symptoms. 

Our study confirmed that increased AWT is associated 
with lower FEV^predicted. This lower FEV^predicted 
depends on the airway size in which the measurements 
of AWT are being performed [12,17] and on the charac- 
teristics of the study population [7,20]. Our population 
consisted of rather healthy elderly males from a ran- 
domly recruited Dutch heavy smoking population and 
still we were able to find significant associations between 
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thicker airway walls and more severe airflow limitation. 
In contrast with the study of Nakano we found a signifi- 
cant negative association between AWT and FEVx/FVC 
illustrating the sensitivity of our method [9] . 

The significant but weak negative association between 
airway wall thickness and emphysema has also been 
reported in other studies [7,13,30] but was not found in 
the study by Nakano [9]. Loss of elastic recoil may con- 
tribute to collapse of the airways resulting in a more 
proximal localization of airways with 3.5 mm internal 
lumen diameter. As these more proximal airways have 
thicker airway walls this phenomenon contributes to the 
weak negative association. Another possible explanation 
for this negative association may be that there are sub- 
jects with predominantly airway wall thickening and 
others with predominantly emphysematous changes. 
Particularly subjects with relatively more airway wall 
thickening are responsible for the negative associa- 
tion and subjects with predominantly emphysematous 
changes do hardly contribute. Apparently, in our popula- 
tion of subjects with normal lung function and with mild 
airflow limitation, the bronchitic phenotype is already 
present in the very early stages of smoking-induced lung 
disease. Discrepancies between the study of Nakano [9] 
and our study may be due to the composition and size 
of the study populations, respectively predominantly 
emphysema versus predominantly healthy smokers with 
respiratory symptoms. 

Importantly, we observed that the contribution of 
AWT 3 5 to airflow limitation was larger than the emphy- 
sema component. Moreover, AWT 35 and emphysema 
together only explained about 40% of the variance in 
FEV!%predicted in this smoking male population. This 
unexpected low contribution of AWT 3 5 and emphysema 
to FEV!%predicted may be that the CT images were 
obtained at full inspiration, while FEV! reflects expira- 
tory airflow limitation. One explanation for this observa- 
tion is that airflow limitation is not only due to reduced 
airway diameter at one level but should be evaluated as 
an integral of all airways at all lumen diameters. This is 
difficult to achieve and therefore we took the smallest 
measurable lumen diameter that provides the largest 
contribution to airflow limitation. A more obvious 
physiological explanation may be the presence of the 
heterogeneity in airway ventilation interrupting the sym- 
metry in parallel airways leading to large clusters of 
poorly ventilated lung units [31]. 

In the univariate analysis, increased AWT 3 5 was asso- 
ciated with respiratory symptoms. However, AWT 3 5 was 
not associated with the presence of any respiratory symp- 
tom in the multivariate analysis after adjustment for 
FEV^predicted, emphysema and smoking behavior. This 
finding corresponds with other studies [7] and is consist- 
ent with the idea that inflammation and airway 



remodeling, associated with chronic bronchitis, is located 
in the more central airways [32]. The study of Martinez 
et al. showed a positive association between airway dimen- 
sions and questionnaires, the BODE index [33] and the St. 
Georges respiratory questionnaire [34] including ques- 
tions about BMI, respiratory symptoms, exercise capacity 
and lung function. Also Camiciottoli et al. found a positive 
association between BODE and airway wall thickness [35] . 
Our study also showed that including respiratory symp- 
toms in the multivariate model with AWT 3 5 and emphy- 
sema has no impact on airflow limitation. 

Lung volume depends on height, weight and sex and 
as a consequence each person has different airway 
dimensions. Therefore, airway dimensions should be 
corrected for lung volume. Actually, volume-corrected 
AWT is the best parameter to use. In this study lung 
volume does not change the multivariate model because 
FEVl%predicted is already corrected for lung volume by 
correcting for patient height. 

It has been shown that the automated imaging pro- 
gram (MEVIS Airway Examiner) based on a method by 
Weinheimer et al. performed much better than the often 
used "full-width-at-half-maximum" method in a silicon 
tube phantom, regarding the blurring effect of CT 
[25,36] . Usually it is better to use sharper kernels for air- 
way quantification. However, it was shown in the study 
by Schmidt et al. that the MEVIS airway examiner pro- 
vides reproducible quantitative results across different 
reconstruction kernels (B30f and B50f) and repeated 
acquisitions [2]. Moreover, the "full-width-at-half-max- 
imum" technique systematically overestimates AWT, 
particularly in small airways [36]. Because low-dose CT 
and the automated imaging program (MEVIS Airway 
Examiner) had not been used previously in clinical stu- 
dies, we firstly optimized our post processing measure- 
ments in smaller cohorts of patients before applying it in 
the population study. We demonstrated that the highest 
number of AWT measurements could be performed on 
airways with an internal diameter of 3.5 mm, irrespective 
of the presence of airflow limitation or emphysema. In 
addition we demonstrated that differences in AWT 35 
are not explained by differences in airway generations. 
Finally, we demonstrated that low-dose CT imaging pro- 
vided similar AWT results as high-dose CT imaging. 

In a non-biased way we were able to evaluate 230 
cumulatively assessed perimeter fractions (APF) per CT, 
ranging from 27-641 APF. In contrast to the commonly 
used pi 10 method, in which a secondary derived variable 
from few, mostly 6 selected airways is used to estimate 
the airway wall thickness [5,6,13-15], we obtained many 
direct airway wall measurements. To our opinion direct 
measurements assessed over all lobes provide a better 
overall reflection of AWT than a limited number of se- 
condary AWT measurements. 
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A limitation of this study is, inherent to general 
population-based studies, that only male smokers with 
mild emphysema and/or airflow limitation were included. 
Strengths of our study is the non-biased way of analyzing 
a high number of airway sections in all lobes that makes 
our approach more suitable for combined airway wall 
thickness and emphysema measurements on one low-dose 
CT scan. Such approach allows monitoring of intervention 
effects on both parameters. This is important when new 
treatment modalities will become available for clinical 
testing. 

In the future further developments may involve mea- 
surements of thickness of airway walls along the full 
length of the bronchial tree at in- and expiration scans. 
Possibly more airflow variability will be explained. 

In conclusion, post processing standardization of large 
numbers of airway wall measurements in all lung lobes 
is feasible, reliable and an useful method to assess airway 
wall thickness. We have demonstrated that increased air- 
way wall thickness contributes more to airflow limitation 
than emphysema in a smoking male population even 
after adjustment for smoking behavior. 

Additional file 



Additional file 1: Table SI. Median AV\^ 3 . 5 and cumulatively assessed 
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FEV^/o predicted. LAA%-950 HU was used to quantify the extent of 
emphysema in this analysis. 
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